The cluster glass state in the two-dimensional extended t-J model 
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The recent observation of an electronic cluster glass state composed of random domains with 
unidirectional modulation of charge density and/or spin density on Bi2Sr2CaCu20s+s reinvigo- 
rates the debate of existence of competing interactions and their importance in high temperature 
superconductivity. By using a variational approach, here we show that the presence of the cluster 
glass state is actually an inherent nature of the model based on the antiferromagnetic interaction 
(J) only, i.e. the well known t — J model. There is no need yet to introduce a competing interaction 
to understand the existence of the cluster glass state. The long-range pairing correlation is not 
much influenced by the disorder in the glass state which also has nodes and linear density of states. 
In the antinodal region, the spectral weight is almost completely suppressed. The modulation also 
produces subgap structures inside the "coherent" peaks of the local density of states. 

PACS numbers: 71.10.-w, 71.27. +a, 74.20.-z, 74.72.-h 
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I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tors (HTS) two decades ago, many anomalous proper- 
ties have been reported. One of the most interesting 
properties is the possible existence of the stripe state 
consisting of one dimensional charge-density modulation 
coupled with spin orderingi^.. The first direct experi- 
mental evidence 4 - of this stripe state as the ground state 
is for the non-superconducting La\ ^N d^^SrQ i 2 CuO^ 
with the hole density about 1/8 per unit cell. Since then 
more evidences of the presence of these have been re- 
ported in other cuprate samples^. However it seems 
that the stripe is much more prominent in the LaSr- 
CuO (LSCO) family near 1/8 doping 8 ' 9 and in partic- 
ular for Lai.s7 5 Ba .i25CuO4 (LBCO — l/8)£ where the 
charge density or spin density modulation could be con- 
sidered as an long-ranged order parameter of a phase 
with broken symmetries. Very recently the high resolu- 
tion scanning tunneling microscopy (STM) has observed 
unidirectional domains with periodic density of states 
modulation in two families of Ca 2 - x Na x Cu0 2 Cl 2 and 
Bi 2 Sr 2 CaCu 2 0$ +5 (BSCCO) 1 ^. These bond-centered 
electronic patterns with a width of four lattice constants 
form the so called electronic cluster glass state with short- 
ranged modulations. Although the modulation is weak, 
it is still surprising to have the rather high supercon- 
ducting (SC) transition temperature in this glass state. 
STM spectra also provided the new puzzling result that 
there are two different types of spectral gaps in those 
systems. One is a larger gap with a broader distribu- 
tion and it seems to be related to the pseudogap. Inside 
this gap, there is a smaller gap or the so called sub-gap 
kink with clear d-wave like spectr a 12 ' 13 . These results 
and the most recent conflicting results reported by angle- 
resolved-photo-emission spectra (ARPES) 14 ' 15 raise the 
question whether there are two energy scales related with 
two different underlying mechanisms separately respon- 
sible for the larger pseudogap and the lower-temperature 
SC transition. But one feature that most experimental 
results agree is the presence of the linear density of state 



(DOS) near the nodal point. It is interesting to note that 
this has also been reported in the non-superconducting 
stripe state of LBCO - l/8±£. 

The observations of these strongly- modulated inhomo- 
geneous states like LBCO — 1/8 with almost zero SC 
transition temperature (T c ) or weakly inhomogeneous 
cluster glass states in BSCCO with quite high T c have 
fueled the idea about the presence of competing inter- 
actions and underdoped HTS is near the boundary of 
two distinct phases with different order parameters. The 
fluctuations 3 of the order parameters in these adjacent 
phases could be the mechanism of high temperature su- 
perconductivity. Thus it seems like a deja vu that after 
two decades we are still faced with a daunting question 
about the appropriate fundamental interactions to model 
and to understand the HTS. Many people are still not 
convinced by all those successe o 17 ' 18 ' 19 claimed by study- 
ing the strong coupling Hubbard model or its equivalent 
t — J model. In this model the nearest-neighbor anti- 
ferromagnetic (AF) spin coupling J is the only relevant 
interaction besides the usual kinetic energy of electrons 
(represented by the t term). There are no other compet- 
ing interactions and certainly no phase boundaries be- 
tween overdoped and underdoped regimes with different 
order parameters to be worried. 

In this paper we will show that interactions repre- 
sented by t and J are actually competing with each other. 
This competition, greatly enhanced by the strong correla- 
tion between electrons, has many spatially heterogeneous 
states almost exactly the same energy as the uniform 
ground state. The ground state could easily tolerate lo- 
cal spatial modulation of charge density, spin density and 
even pairing amplitude without much an effect on its SC 
order parameter. The presence of these very low-energy 
cluster glass states with a random pattern of short-ranged 
modulation is an inherent nature of the t — J model. The 
selection of a particular local electronic pattern as ob- 
served in scanning tunneling spectroscopy (STS) for HTS 
is most likely determined by the effects of impurities, de- 
fects and electron-lattice interaction, etc. Random dis- 
tribution of impurities and defects will not produce long- 
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ranged modulation in the sense of charge-density-wave 
order or spin-density-wave order unless these is a very 
strong lattice distortion as demonstrated^ by the struc- 
tural transition observed in LBCO — 1/8. Hence our re- 
sult for the extended t — J model shows that for weakly 
inhomogeneous cuprates like BSCCO there is no need to 
introduce other strong competing interactions with new 
broken symmetry phases and long-ranged order parame- 
ters to produce the observed glass states. The extended 
t — J model is adequate to explain many of the experi- 
mental observations. Depending on the particular type of 
modulations, the local DOS could either have a node with 
linear DOS or without a node. These local modulations 
also produce the sub-gap structures. In addition, these 
cluster glass states have almost completely suppressed 
the quasi-particle spectral weight near the antinodal re- 
gion as observed by ARPES for BSCCO compounds 20 . 

Before we start with the discussion about the clus- 
ter glass states, we should point out that the competi- 
tion between the kinetic energy represented by t and the 
magnetic energy represented by J is not a new idea. It 
is the main reason for the uniform-state phase diagram 
obtained by the theory of the resonating-valence-bond 
(RVB) stat o 17 ' 18 i 19 . The J term prefers the formation of 
spin pairing, the so called RVB singlet, or the long-ranged 
antifcrromagnetism. Because of the strong correlation a 
spin can only hop by exchanging position with a hole, 
hence the kinetic energy is proportional to the density of 
holes. But the more hole the system has, the less number 
of spin there is. Consequently the pairing is also reduced. 
Thus the pairing amplitude reduces when hole doping in- 
creases. Magnetic energy is reduced while kinetic-energy 
gain increases. It is shown below this competition also 
happens in individual unit cell. 

On the theoretical side the stripe state also has a long 
history. It was first founded almost two decades ago 
in the mean-field treatment of the Hubbard model 12 
although the validity of the method in treating strong 
Hubbard interaction is questionable. Then the extended 
t — J models were studied by several numerical methods 
such as the exact diagonalization method 2 -^, the density 
matrix renormalization group metho d 22 ' 23 and the varia- 
tional Monte Carlo (VMC) method^!. But the results 
are inconsistent. There are indications that stripe is un- 
stable when the second neighbor hopping, t', is included 
in the t — J model. However, a later VMC study of the 
t — t' — J model 2 ^ indicates that the stripe has about 
1% lower energy than the uniform RVB based e?-wave SC 
state for most of the negative values of t' except when 
— t'/t is less than 0.1. Similar results were also reported 
for the Hubbard model^. These results pose a new direct 
contradiction with the experimental findings. Many ex- 
perimental and theoretical studies have found that —t'/t 
is smallest for LSCO family 2 ^ 28 ' 29 , thus the stripe state 
should not be favorable. Yet as mentioned above, among 
all the cuprates LSCO family has the most solid evidences 
for the existence of stripe. There are also other prob- 
lems with the proposed stripe states, such as suppressed 



pairing correlation and absence of the V-shape DOS at 
low energy, and we will discuss them below. All of the 
previous works concentrated on studying periodic stripes 
with a long-ranged order, it is unclear if the result will 
hold for the cluster glass state. There are also differ- 
ent kinds of stripes to be considered. It is possible to 
have only charge density modulation, or only spin den- 
sity modulation, or only pairing amplitude modulation 
or with the linear combination of any two or all three of 
the m 30 ' 31 . Furthermore the relation between these mod- 
ulations could be correlated or anti-correlated. All these 
issues are addressed below and their results are compared 
with the experiments. 

II. THE STRIPE-LIKE STATES BY THE 
VARIATIONAL MONTE CARLO METHOD 

We consider the extended t — J Hamiltonian, 

i,j,<r <i,j> 

which has been used to describe many physical systems 
in the high-temperature superconductors^. The hop- 
ping amplitude Uj =t,t', and t" for sites i and j being 
the nearest-, the second-nearest, and the third-nearest- 
neighbors, respectively. We restrict the electron creation 
operators c\ a to the subspace with no-doubly-occupied 
sites. Si is the spin operator at site i and < i,j > 
means that the interaction occurs only for the nearest- 
neighboring sites. In the following, we mainly focus on 
the case t" = — 1'/2 and J/t — 0.3 at hole doping 1/8. 

We shall follow the work by Himeda et al£5. to con- 
struct the variational wave functions. In the mean-field 
theory we assume a local AF order parameter, the stag- 
gered magnetization mj , and nearest neighbor pairing or- 
der parameter Aij. Thus the effective mean- field Hamil- 
tonian is reduced to 

where the matrix elements 

H ijrT = - I it, E +t '» E +t " E ^'•*'+/ 3 

\ p=N I3=NN P=NNN/ 

+ ( Pi - Mw + a(-l) Xi+yi mi ) 6 jti , (3) 
D v = A *A*+/3- ( 4 ) 

/3=N 

Here j3 — N, NN, and NNN correspond to the nearest- 
, the next-nearest, and the third-nearest-neighbors, re- 
spectively, and a =T(1) or |(-1). The local charge den- 
sity is controlled by pi and [x v is the variational param- 
eter for the chemical potential. For periodic stripes we 
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assume charge density pi and staggered magnetization 
rrii are anti-correlated, i.e. there are more holes at sites 
with minimum staggered magnetization. For simplicity, 
we assume these spatially varying functions with simple 
forms: 



Pi = p v cos[47r<5 • (m - yo)], 



rrii = m v sin[27r6 • - y )], 



(5) 



(6) 



where <5 is the doping density and is 1/8 in this paper. 
We have also taken the lattice constant to be our unit. 
Here we assume the stripe extends uniformly along the x 
direction, yo = (1/2) corresponds to the site- (bond-) 
centered stripe. In this paper, we will only focus on the 
bond-centered stripe since the variational energy differ- 
ence between the site- and bond-centered stripes is very 
small (not shown) within the finite cluster— 

Besides staggered magnetization and charge density, 
the nearest neighbor pairing Ay can also have a spa- 
tial modulation. There are several different stripes we 
can choose. If Ay has the same period 1/5 as the stag- 
gered magnetization but it is n phase shifted, this is 
the so called "antiphase" stripe studied by a number of 
group o 25 i 33 ' 34 i 35 . In this stripe state, the bond-average 
Ay is zero and there is no net pairing. Hole density is 
maximum at the sites with maximum pairing amplitude 
|Ajj | and minimum staggered magnetization \rrii\. 

Another more general choice is to have the spatial vari- 
ation of A,-, of the form, 



A 



Af cos[4^ • (y t - y Q )] - A 



c 



A M+fi = -A* 1 cos[47r5 • ( w - yo ) + 2ttS] + A^ . (7) 

The bond-average Ay is determined by the constant A„ . 
If both A^) f and A^ are positive, then the hole density is 
maximum at sites with smallest pairing amplitude |Ay| 
and smallest magnetization |m<|. This is similar to the 
phase diagra m 36 ' 37 predicted by the uniform RVB and 
AF states, when hole density is small both staggered 
magnetization and pairing amplitude are larger. Thus we 
will denote this state as the AF-RVB stripe. For AF-RVB 
stripe the period of Ay is the same as the charge-density 
modulation pi. This period, 1/2S, is only half of the 
period for the antiphase or tt phase stripe. Besides the 
antiphase stripe and AF-RVB stripes, we could also have 
the AF stripe without both A^f and A^ or the charge- 
density stripe without any staggered magnetization but 
with pairing amplitude modulation. 

In general we have total seven variational parameters 
p v , t' v , t", p v , m v , A^f , and A^ with t v set to be 1. Once 
these parameters are given, we diagonalize the mean-field 
Hamiltonian in equation @. By solving the Bogoliubov 
de Gennes (BdG) equations 



and then obtain N positive eigenvalues E n (n = 1 — 
N) and N negative eigenvalues E n with corresponding 
eigenvectors (u™,u™) and (u",vf). The eigenvectors are 
used to construct the mean-field wave function^ by 
using Bogoliubov transformation 
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The trial wave function P\ip) with a Gutzwiller projec- 
tor P can be constructed by creating all negative energy 
states and annihilating all positive energy states on a 
vacuum of electrons |0). Then we formulate the wave 
function in the Hilbert space with the fixed number of 
electrons N e , 



|$) =PP N M = PPN e Hlnll\0) 



(10) 



A/g/2 
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where [7y = and Vij — v*. We optimize the varia- 
tional energy E = ($|iJ|$)/($|$) by using the stochas- 
tic reconfiguration algorithm 38 . Additionally, to reduce 
the boundary-condition effect in numerical studies^, we 
average the energies over the four different boundary con- 
ditions which is periodic or antiperiodic in either x or y 
direction. 

FigfT] shows the t' /t dependence of the variational en- 
ergies for the hole density 5 = 1/8. Firstly, we have 
shown that the AF-RVB stripe state (the empty trian- 
gles) becomes more stable than uniform <i-wave RVB 
state (the filled triangles) as decreasing t'/t further from 
—0.05. It is worth noting that the AF-RVB stripe state 
for t'/t < —0.05 has the vanishing pairing parameters 
A^ and A^f but large m v and finite p v . Therefore, we 
can consider this stripe state without SC order as the an- 
tiphase AF stripe state 1 . The hole density and the stag- 
gered magnetization are plotted as a function of positions 
for a typical AF stripe (filled circles) in FigJ^a) and (b) , 
respectively. The AF stripe pattern with a very large 
hole-density variation is essentially a nano-scale phase 
separation with hole-rich and hole-poor regions alternat- 
ing. This is consistent with what has been obtained by 
Himeda et al^- although they did not include t" which 
is set to be —t'/2 here. 

In our previous calculations to study the possibility 
of phase separation in the t — J model^ 6 -, we found the 
tendency to overestimate the strength of the pairing of 
holes. If we reduce this strength by going beyond the 
simple trial wave functions used in our discussion above 
we could push the phase separation boundary to a much 
higher value of J/t. A simpler way to make this ad- 
justment is to introduce the hole-hole repulsion Jastrow 
factoiMdi: 
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FIG. 1: Variational energies of the bond-centered AF-RVB 
stripe state and uniform d-wave RVB state as function of t'/t. 
The data have been optimized for 1/8 doping on the 16 x 16 
lattice system. The filled (empty) symbols represent the ener- 
gies of uniform d-wave RVB state (AF-RVB stripe state). The 
triangles (circles) indicate the wave functions without (with) 
the hole- hole correlation of equation (|11[) . The asterisks cor- 
respond to the random stripe states with the Jastrow factor of 
equation (|11[) . Inset: The hole density of the random stripe 
state optimized in the case of (t',t" , J)/t = (—0.3,0.15, 0.3). 



with 



r " = \/ snr I j~ u '- -''.'O + sin2 [j^{yi-yj))i 



where n^l = 1 — Y< a cl a Ci a . The three parameters vp of 
(3 = N, NN, and NNN are for short-ranged hole-hole 
repulsion (attraction) if these values are less (greater) 
than 1. The factor is for long-ranged correlations^, 
and it is repulsive if a is positive. L x and L y are the 
number of sites in the x and y direction, respectively. 
We have found that the variational energy of the peri- 
odic stripe state without including the Jastrow factors is 
very sensitive to boundary conditions. However the Jas- 
trow factor is capable of reducing the dependence of the 
boundary conditions. 

As shown in Fig[TJ when the hole-hole repulsion of 
equation Jll]) is included in our VMC calculation, both 
the uniform d-wave RVB state (filled circles) and the AF- 
RVB stripe state (empty circles) have gained significant 
energies. In particular, the uniform d-wave RVB state 
has gained about 2 ~ 3% in total energy. It indicates that 
for the extended t — J model the holes are not bounded as 
tightly as estimated by the simple RVB trial wave func- 
tions. It also implies that the nano-scale phase separation 
of the AF stripe will not be preferable. Indeed, the low- 
est energy stripe state does not favor the AF stripe but 
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FIG. 2: The profiles of (a) hole density, (b) staggered mag- 
netization and (c) pair-pair correlation function with R y = 8 
for the optimized states at 1/8 doping in the extended t — J— 
model with t'/t = -0.3 and J/t = 0.3. The filled (empty) 
circles correspond to the AF-RVB stripe state without (with) 
hole-hole repulsion. In (c), the dashed (dashed-dotted) line 
indicates the uniform d-wave RVB state without (with) hole- 
hole repulsion. The red solid line corresponds to the random 
stripe state with hole-hole repulsive Jastrow correlation. All 
quantities are calculated on a 16 x 16 lattice system with the 
periodic and antiperiodic boundary conditions along the x 
and y directions, respectively. 



the AF-RVB stripe with a much more reduced staggered 
magnetization as shown in FigfJUb) and a larger pairing 
parameter A„ . Now the hole density has a much smaller 
variation as shown by the empty circles in Fig[2ja). 

It is surprising to find out that for all t'/t, the AF- 
RVB stripe state is almost degenerate in variational en- 
ergy with the uniform d-wave RVB state. This is quite 
remarkable as the two trial wave functions are very dif- 
ferent. As an illustration we show in FigfJJa) and (b) 
the variation of the hole density and the staggered mag- 
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netization along the modulation direction, respectively, 
for an AF-RVB stripe state (empty circles). It should 
be noticed that not only the hole density is completely 
uniform for <i-RVB state but there is also no staggered 
magnetization at all. Instead of periodic stripes, we have 
also examined the stripe state with 4x4 patches in the 
16 x 16 lattice system. For each patch, we choose a di- 
rection of the stripe, x or y, randomly. We consider this 
state as random stripe state. For simplicity, we still use 
the same equations (j5 l7lll[) . The hole-density modula- 
tion of a random stripe state has been shown in the inset 
of FigfT] As shown in Fig[TJ the random stripe states with 
the Jastrow factors of equation (fTTj) also have the opti- 
mized energies almost identical to uniform d-wave RVB 
and AF-RVB stripe states even though we have not op- 
timized parameters on very bond or site independently. 
These optimized random stripe states have finite m v and 
A^ but smaller A* f (less than one third of A%). 

We believe this energy degeneracy is caused mainly by 
two reasons. The first reason is that the terms in the 
t — J Hamiltonian are all local within nearest neighbors 
or next nearest neighbors. The second reason is that the 
hopping terms and spin interaction not only are of the 
same order of magnitude but they are competing against 
each other. We have found the energy competition be- 
tween the kinetic energy and spin interaction is very ro- 
bust among different states (not shown). Their competi- 
tion is enhanced by the no-doubly occupancy constraint 
as the presence of holes will suppress the spin interaction 
to zero. Thus it is possible to have locally different spin- 
hole configurations with different emphasis on the kinetic 
energy or the spin energy. Some of these patterns have 
lower kinetic energy but higher magnetic energy than the 
uniform tf-RVB state and some with opposite energetics. 

There are at least two important implications of this 
energy degeneracy. The first one is that the inhomoge- 
neous states are quite robust in the extended t—J models 
without the need for introducing any other large inter- 
actions. States with different local arrangement of spin 
and holes may have very similar energies. The second 
implication is that any small additional interaction could 
break the degeneracy. If we consider the " realistic" situ- 
ation of cuprates with large numbers of impurities, disor- 
ders and electron-lattice interactions, the inhomogeneous 
states could be much more numerous and complex than 
we have expected. Materials made with different pro- 
cessing conditions could also produce different inhomo- 
geneous states. Since all these are presumably secondary 
interactions smaller than the dominant t and J, the mod- 
ulations of charge density, staggered magnetization and 
pairing are expected to be small. If one of the inter- 
actions, like electron-lattice interaction, becomes quite 
strong as observed in LBCO — 1/8^, then the modula- 
tion will also become larger and longer ranged. 

We have also investigated the pair-pair correlation 
function for the optimized states with/ without hole- hole 
repulsive Jastrow factors in the case of (t',t", J)/t — 
(—0.3,0.15,0.3). The singlet pair-pair correlation func- 



tion along the modulated direction (y direction) is de- 
fined as 



1 ^|($|At(r)A x (r + J R y )|$) 



, (12) 



c ii c l+xi creates 



where r = (x, y) and Aj.(r 

a singlet pair of electrons among the nearest neighbors 
along x direction for each site r. Here, we focus on the 
long range correlation, and thus set R y — 8 to be the 
largest distance on 16 x 16 lattice system. In FigJ^c), it 
is shown that the hole-hole repulsion suppresses the long- 
range pair-pair correlation about three times in the uni- 
form cZ- wave RVB state. We have also found that without 
hole-hole repulsive correlation, the long-range pair-pair 
correlation is much more reduced in the AF-RVB stripe 
state than uniform d-w&ve RVB state, because the AF- 
RVB stripe state has almost vanishing A^ . Due to large 
m v , this AF-RVB stripe state shows much larger ampli- 
tude of the P xx modulation. However, after considering 
the hole-hole repulsion, the AF-RVB stripe and uniform 
d-wave RVB states have almost the similar magnitude 
of pair-pair correlation as shown by the empty circles 
and dashed-dotted line in FigfJJc), respectively. For the 
random stripe state shown in the inset of FigfT] the aver- 
age value P xx for the whole system is shown as the red 
solid line in Fig^c). It is essentially the same as the 
value of uniform d-wave RVB state and the periodic AF- 
RVB stripe state. Although the staggered magnetization 
for the random stripe state has larger variation than the 
periodic stripe state shown in FigJ^b), the pairing cor- 
relation is unchanged. These results indicate that the 
long-range pair-pair correlation is mostly determined by 
the value of A^ and is rather insensitive to the modu- 
lation of the hole density and staggered magnetization. 
Thus we also expect to have a robust cZ-wave node ob- 
served by recent experiment s 10 ! 13 . This will be shown 
below. 



III. DENSITY OF STATES BY THE 
GUTZWILLER APPROXIMATION 

According to the VMC calculation for the extended 
t—J model, it is likely that there are a number of inho- 
mogeneous states close in energy to the uniform ground 
state. Then, some sort of small perturbation may choose 
a particular stripe state as the ground state. Assum- 
ing such a situation, here we regard a stripe state as the 
ground state, and consider the projected quasi-particle 
excitation spectra. However, calculation of the excited 
states by the VMC method is computationally very ex- 
pensive, and it takes too much time to investigate wide 
parameter range to obtain general properties of the stripe 
states. Furthermore, one can take only a limited system 
size and it is difficult to obtain dense spectra. 

Therefore, as a first step, we use a Gutzwiller mean- 
field approximation for this purpose. The minimization 



6 




FIG. 3: Spatially averaged local DOS for the random stripe 
states calculated by the non-self-consistent BdG equation, (a) 
For A? / 0, we use parameters optimized by the VMC, t v = 
-0.35, i" = 0.16, m = 0.15, p = 0.03, = 0.28 , Af = 
0.02, but /i = — 0.875£„ is adjusted to realize 1/8 filling, in 
units of t v . (b) The same parameters except for Ajf = 0. 



of the total energy yields a BdG equation 9 . Usually 
the parameters in the BdG equations are solved self- 
consistently to find an optimal solution. However, since 
we already have the assumed inhomogeneous ground 
state here, we shall use parameter sets obtained from 
VMC results and diagonalize the BdG Hamiltonian only 
once, instead of solving self-consistently. Furthermore, 
for convenience, we slightly modify the Gutzwiller pro- 
jection by attaching fugacity factors. Namely, we assume 
that P\\ip) is the ground state and that P\j^\ip) are the 
excited states, where P\ = P ]J icr A™j" , and \ ia is a fu- 
gacity factor to impose the local electron density conser- 
vation for each spin, namely, W ^p^) W = for 

any i and a. The quasi-particle operators 7^ are obtained 
by solving the BdG Hamiltonian. In addition, here we do 
not take into account the Jastrow factor as it only affects 
the hole-hole correlation slightly but not the local DOS 
studied below. Then, under the assumption of the non- 
self-consistency, the BdG Hamiltonian is represented by 
equation (J2J). With this formulation, P\jl]ip) of differ- 
ent n are approximately orthogonal to each other—, and 
thus we expect that it is suitable to use P\ instead of P 
for our purpose here. Since the result presented below 
is qualitatively not very sensitive to small change of the 
parameters, we expect that such a modification of the 
projection should not affect the results qualitatively. 

Then, by taking the most dominant terms, the local 
DOS is represented by 

N<t(R,u>) = fa J2 Kl 2(5 ( w ( 13 ) 

n 

N l (R^)=g t Rl J2K\ 2 ^ + E n ), (14) 



where index n runs for both positive and negative eigen- 
values. Note that only the position dependent constant 
g Ra = (1 — riijf — Uri)/(1 — TiRa) is multiplied in front 
of the local DOS by the standard BdG formalism. Since 
the result of the site-centered stripe is very similar to 
that of the bond-centered stripe, we show only the lat- 
ter. Here we shall only discuss our results for the ran- 
dom stripe state. For the random stripe state consid- 
ered above, each randomly oriented domain is assumed 
to have the same parameters so that the VMC calcula- 
tion is possible. Ideally we should have optimized these 
variational parameters on every bond or every site. Then 
we expect to have a much broader distribution of these 
parameters. To simulate this effect, we simply replace 
each Aij by (1 + £^ 3 -)Ay, where £y is a random vari- 
able which has the Gaussian distribution around with 
the standard deviation of 1. We use a supercell of size 
32x32 sites, and the same configuration is repeated as 
20x20 supercells to obtain the local DOS. The Fourier 
transform with respect to the supercell index is similar to 
a system of small clusters with many twisted boundary 
conditions. 

As shown in FigJ3Ja), the spatially averaged 
J2 a N a (R,u)) for the AF-RVB stripe (A£ ^ 0) has a V- 
shape at low energy. On the other hand for an antiphase 
stripe which has = there is no V-shape as shown in 
FigEtb). The energy level is broadened with a width of 
0.05^. In FigEJb) , there is a very small dip at E = 0, it 
would be bigger if A* f increases. Our results are consis- 
tent with the very recent report by Baruch and Orgad 31 . 
In general, there is no V-shape DOS for the antiphase 
stripe. For the same system as FigU^a), in Figdjwe plot 
the position dependence of local DOS at randomly cho- 
sen 8 sites. The low energy spectra seem less influenced 
by the disorder than high energy. This result shows that 
the node and the low energy V-shape DOS are robust 
against this kind of inhomogeneity. This is possible be- 
cause nodal fc-points do not have many states to mix with 
and also the suppression of impurity scattering 50 . The 
sub-gap structure o 12 ' 13 also seem to be quite apparent. 
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If we switch off the random variables the gap vari- 
ations from site to site are small. These gap variations 
grow larger as distributions of £y become wider. 

The key to understand the absence of V-shape in the 
Local DOS for the antiphase stripe is the Fourier trans- 
form of the modulated Ay term written as, 

A^5>sfc x (e- ie 4 +giT ct_ M + e*4_ g , T c t _ Jfe)i +h.c.) 

k 

_ Af 
2 



k L 



e -i{k v +q y ) \ R i(- 



U t 



-k-q,^-k,[ 



h.c. 



(15) 



where 9 — for site-centered stripes and 9 = q y /2 for 
bond-centered stripes; q — (0,7r/4) for the antiphase 
stripe, and q = (0, tt/2) for the AF-RVB stripe. What 
is important here is that it contains only pairing with 
nonzero center-of-mass momentum as the FFLO stated. 
In the case of zero-momentum pairing as the conven- 
tional BCS theory, the spin-up electron band couples 
with the spin-down hole band (the upside-down down 
electron band) . These bands intersect at the Fermi level, 
and a gap opens if Afe 0. In the case of finite-q pair- 
ing, however, the spin-up electron band couples with ±q 
shifted spin-down hole bands, and thus the band inter- 
sections occur not at the Fermi level. Therefore, a gap 
does not open at the Fermi level. The constant A^ 7 term 
which forms the usual Cooper pair is necessary for having 
the node and the V-shape DOS. 

To compare with ARPES experiments, A(k,u>) is also 
calculated. Since A(k,u>) is regarded as the local DOS in 
the A;-space, Let us take the Fourier transform of renor- 
malized u R , v R , namely, 



i 



(16) 



Then, A a (k,Lu) is written as 



A t (fc,,) = ^K| 2 %-£„), (17) 

n 

4 1 (fc, w ) = ^K| 2 % + £ I1 ). (18) 

n 

Fig [5] shows A a (k, uj) of the random stripe state 
with A^ 0, where each of the ^-function spectra are re- 
placed with a Lorentzian distribution with 5E = 0.05^. 
In comparison with the uniform d-wave RVB state, the 
spectral weight around the antinodal region is much 
strongly reduced than near the node. 

IV. CONCLUSIONS 

In summary, we have used a variational approach to 
examine the possibility of having inhomogeneous ground 
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FIG. 5: 57J<r A a (k, E) for the random stripe state with 
Ac 7^ 0. Dotted lines are for uniform d-wave RVB state 
{p v — 0, m v = 0, A* 7 = 0). (a) Nodal and (b) near antinodal 
region. 



states within the extended t — J model with 1/8 doping. 
We considered states with spatial modulation of charge 
density, staggered magnetization and pairing amplitude. 
Besides the antiphase or inphase stripes considered by 
many groups we have proposed a new AF-RVB stripe. In 
this stripe state, we assume there is a constant pairing 
amplitude besides the various modulation. In addition 
to considering states with periodic stripes we also con- 
sider random stripe states to simulate the cluster glass 
state observed by experiments 10 . By improving the trial 
wave functions with the introduction of hole-hole repul- 
sive correlation, we have greatly improved the variational 
energies by several percents for both uniform RVB d-wave 
SC state and states with periodic AF-RVB stripes for re- 
alistic values of i'/t. Most surprisingly the random stripe 
state essentially also has the same energy as the uniform 
state in spite of our oversimplified assumption that all 
the stripe domain has the same patterns of modulation 
instead of each site or bond with different values. This 
random stripe state also has about the same long-range 
pair-pair correlation as the uniform or periodic stripe 
state even though there are significant staggered mag- 
netization and charge variation from site to site. Then 
we also examined the local DOS and the spectral weight 
of the random stripe state by using Gutzwiller approx- 
imation. We found the V-shape DOS and the node are 
still present at every site. The local DOS measured at dif- 
ferent positions shows a broad variation of the gaps and 
also it has sub-gap structures seen in experiment o 12 ' 13 . 
The spectral weight at the antinodal direction is negli- 
gibly small but finite around the node. All these results 
are quite consistent with experiments reported for cluster 
glass state in BSCCOi£. 

Our result also resolves an inconsistency with ex- 
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periments derived from previous theoretical calculations 
without including the hole-hole repulsion in the trial wave 
functions. Stripe is neither stabilized nor destabilized by 
the long range hopping. In fact, due to the competition 
between the kinetic energy gain and magnetic interac- 
tion, it is very natural to have the spatial modulation, in 
periodic or random configuration, of charge density, mag- 
netization and even pairing amplitude. The constraint of 
disallowing doubly occupation of electrons at each lattice 
site has significantly enhanced the competition. Many 
local arrangements of spin and hole configurations could 
give almost identical total energy as the uniform solution. 

Recently, Capello et alM- have also found that the en- 
ergy of the periodic RVB stripe state is very close to 
that of the uniform RVB state by using a variational cal- 
culation. They have considered several possibilities for 
the stability of the RVB stripe state, such as lattice dis- 
tortion, t'-effect, and long-range Coulomb repulsion with 
the conclusion that uniform RVB state is still the lowest 
energy state. This is very consistent with our conclusion 
although we have included the antiferromagnetic order in 
the stripe state. The issue about whether AF is present 
in the stripe will be discussed in the future. They have 
not considered the cluster glass state with random AF- 
RVB stripe domains which is also a good candidate for 
the ground state. 

The presence of inhomogeneous or cluster glass states 
is apparently a very natural consequence of the t — J 



model. There is no need for introducing additional in- 
teractions to generate such states. In fact we showed 
that the RVB state with a finite constant pairing is quite 
compatible with the local variations of charge density, 
magnetization and even pairing amplitude. As long as 
this modulation is not overly strong, the superconductiv- 
ity still survives as the node and V-shape DOS are still 
present. In a realistic material, other interactions such as 
impurity, disorder, and electron-lattice interactions, etc., 
no doubt will help to determine the most suitable local 
configuration of spins and holes but they will not produce 
a globally ordered state unless there is a very strong and 
dominant interaction like the electron-lattice interaction 
seen in Lai~ x Ba x CuO\ at 1/8 doping. The verification 
of this is left for future work. 
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